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1. INTRODUCTION 

Chaotic systems are widely used in wireless communication systems due to their non-periodicity, the 
randomness property, and the high sensibility to the initial values [1], [2]. Among all types of chaotic based 
modulation systems, a non-coherent chaotic modulation scheme has been highly interested by researchers due 
to the simplicity of the design and removing the wanted chaotic synchronization circuit at the reception. The 
famous practical circuit of non-coherent chaotic modulation shift keying is the differential chaos shift keying 
(DCSK) system [3], in which the system is designed with a simple circuit with an appropriate bit error rate 
(BER) under multipath fading channel without needing channel state information at the reception [4]. 
Nevertheless, in classical DCSK systems, only half of the bit period is consumed on carrying the data sequence 
[5]-[7], and the multicarrier based DCSK structure has been proposed by Kaddoum et al. [5], such as the 
orthogonal frequency division multiplexing (OFDM), to achieve higher information rate transmission and 
improve the energy efficiency of the DCSK system, which is also named multicarrier-DCSK (MC-DCSK) or 
OFDM-DCSK system. The subcarriers are divided into multiple groups, only one chaotic reference is sent for 
each group [8]. OFDM-DCSK is a non-coherent chaos-based communication structure. Improving data rate 
and saving bit energy are the most important advantages of this system, it’s also simple to design as it does not 
require RF delay, and adding to the mentioned advantages, there is no need to use a chaotic sequence generator 
at the receiver end. The high PAPR of the transmitted signal is the major drawback of OFDM-DCSK. The RF 
power amplifier (RF-PA) is affected by the high PAPR, which gives rise to signal distortion. A high PAPR 
degrades the operation of the analog-to-digital and digital-to-analog converters and minimizes the performance 
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of the RF-PA, this distorts the transmitted OFDM-DCSK signal, which leads to poor BER performance of the 
whole system. 

The PAPR is defined as the ratio of maximum power to the mean power of the transmitted OFDM- 
DCSK signal during one symbol. It’s also expressed in units of dB. PAPR occurs in a multicarrier system when 
the different subcarriers are out of phase with each other. The problem of PAPR was investigated by authors 
and they have suggested PAPR reduction algorithms to improve the performance of the MC-DCSK systems 
[9]. They insert dummy sequences in accordance with the inverse fast fourier transform (IFFT) [10] needs no 
receiver modifications. Also, several techniques were suggested to overcome the issue of PAPR for OFDM 
systems, like palm date leaf clipping [11] and partial transmit sequences (PTS) [12]. A standard block code 
[13] was used for error correction and for OFDM PAPR reduction. In 2011, Lain et al. [14] proposed a PTS 
technique by using real-valued genetic algorithms to reduce the PAPR and computational complexity. PAPR 
reduction methods for both MRT (maximum ratio transmission) OFDM and EGT (equal gain transmission) 
OFDM systems were also suggested in 2014 [15], based on the obtained findings. It's worth noting that the 
suggested method can enhance both PAPR and the bit error rate in MRT OFDM systems, whereas it improves 
PAPR but only marginally decreases the bit error rate in EGT OFDM systems. 

Kang et al. [16] proposed a hybrid approach of WHT precoding and companding, this approach is 
reducing the PAPR without affecting the performance of BER and PSD. Liu etal. [17] presented interferometry 
code-aided to reduce the PAPR of OFDM-DCSK. Each chaotic chip in this system is distributed over all 
possible subcarriers via CI codes, thereby decreasing the coherence among the signals and lowering the PAPR. 
The selective mapping technique (SLM) is one of the most effective techniques that has been used to reduce 
PAPR. Jayalath et al. [18] reduced PAPR by using Newman phase sequences. Goel et al. [19] presented the 
SLM technique to reduce the PAPR by using monomial phase sequences, which also reduced the computational 
complexity. Goel et al. [20] implemented a mixture of M-ary chaotic sequences and a mapping system. The 
major benefit of this technique is that it does not require side information transmission and gives good BER 
performance with good PAPR reduction. they [21] presented a combined PAPR reduction technique by using 
Riemann SLM and DCT. Mobini and Zahabi [22] presented a merit factor SLM technique to improve the 
PAPR reduction as compared with the conventional SLM technique. Different companding techniques were 
investigated, such as hyperbolic, u-Law, and A-Law, which were compared with the palm clipping technique 
to study the effects of these techniques on the PAPR of the OFDM-DCSK system and the BER performance 
[23]. The contribution of this paper is to reduce the PAPR value of OFDM-DCSK system by utilizing SLM 
technique using Hadamard, Riemann, and enhanced matrices as phase sequences and compare it with 
conventional SLM technique. 


2. SELECTIVE MAPPING (SLM) PAPR REDUCTION TECHNIQUES BASED OFDM-DCSK 
The PAPR is defined as the ratio of an OFDM-DCSK symbol's maximum power to the average power 
of that symbol. 


Max (|s(n)|?) 


PAPR = TESE 


,when0 <n < N-1 (1) 


where s(n) is the transmitted signal of the OFDM-DCSK system and N is the length of the sequence. PAPR 
is determined by the input data. PAPR is a random variable since the input data is a random variable. The 
likelihood that instantaneous power is above a given threshold value can be calculated using the amplitude 
distribution of the OFDM-DCSK signal. The CCDF for various PAPR threshold values is used to calculate it. 
The efficiency of the PAPR reduction technique is conventionally used as a performance measure. The chance 
of an OFDM-DCSK block being more than a certain threshold value is defined as: 


CCDF(PAPR,) = Pr(PAPR > PAPR,) (2) 


The threshold value and probability are defined in the above equation as PAPR, and Pr respectively. 

SLM is a PAPR reduction technique that converts an OFDM-DCSK signal into many separate signals 
using phase sequences and transmits the signal with the lowest PAPR. For the OFDM-DCSK system, SLM is 
one of the distortion less PAPR reduction techniques. SLM needs the index of the chosen signal to be 
transmitted along with the OFDM-DCSK symbol. The sum of the change in PAPR output is determined by the 
number of candidates U and the kind of step sequences used. The reduction in PAPR is proportional to U, but 
increasing U raises side knowledge overload, which is expressed as (log2U). When it comes to lowering PAPR, 
step sequence selection is crucial. The block diagram of the OFDM-DCSK selective mapping technique is 
shown in Figure 1. 
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Figure 1. Block diagram of OFDM-DCSK transmitter for selective mapping technique (SLM) 
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In the first, the chaotic sequence is generated using second-order Chebychev Polynomial (CPF), which 
is given X44, = 1 — 2xf so that the | — th reference sequence is defined as X,=[Xj9, Xi) +. ++ a Xg- 
where f is the spreading factor. The mean value for the generated sequence is zero and the mean square value 
is unity, E [X,]=0 and E [X7]=1. The l— th parallel data sequence, bım, m = 1,2,..,M € {—1,1}, is 
generated from the parallel mapped stream bits using a serial to parallel converter, where M is the length of a 
parallel data sequence that are sending in each frame. The parallel data sequence is multiplied by a chaotic 
reference sequence and then IFFT is taken for each k — th vector (k = 0,1, ..., P — 1) after multiplying it by 
the u — th phase vector, p“, u = 1,2,..,U. There are 6 number of IFFT transforms that are required to 
complete one frame. The l — th transmitted OFDM-DCSK signal and the u — th phase code can be expressed 
as: 


j2mmn 


1 1 
si (n,k) = Jr um=1 bimXix * el N ) X Pm.Lk + Ja tuk £ Polk» (3) 
0< n < N-1, k=0,1,..,B-1 


The l — th transmitted OFDM-DCSK signal which has the lowest PAPR value will be selected to be 
transmitted along with the & index and the l — th transmitted signal is & (n, k) . N is the FFT size, N = M + 1. 
We assume here that the cyclic prefix is not considered. The k — th received signal is expressed as [5], [24]: 


n(n, k) = Xia yis (n — Tik) +w(n,k), 0 < n < N-1,k=0,1,..,B-1 (4) 


Where y; and T; are the Rayleigh channel coefficient and the corresponding time delay for the i — th path, 
respectively. L is the number of paths and w, is the 1 — th AWGN with zero mean and variance of N,/2. The 
Rayleigh pdf of y; is written as [5]: 


v2 
fpo) = je 27, v>0 (5) 


Where ø is the standard deviation of the distribution that is greater than zero. At the receiver, non-coherent 
detection is used. There is no channel estimation required as in the OFDM system and no RF delay required as 
in the DCSK system. The l — th received signal, r, is converted to a parallel sequence by using a serial to 
parallel converter. The l — th parallel sequence, ñi, i = 0,..,M, are passed through the FFT transform to 
obtain the m — th and k — th index of the OFDM-DCSK demodulated signal at the l — th frame as: 


—j2mmn 


Ry(m, k) = ENAC K) * ew) spi p,m = 0,1,..,M,k =0,1,..,B —1 (6) 


The first subcarrier (zero frequency) contains the received reference chaotic sequence, R,(0,k) and 
the remaining subcarrier contains the M received information chaotic sequences, Rim, k),m = 1,..,M. 
There are M correlators are produced from these sequences by multiplying the reference chaotic sequence by 
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the m — th information sequence and summing over p period. Then the output of the m — th correlator is 
given by: 


Qim = ERZE R,(0,k) * B,C, k), m = 1,2,...,M (7) 


For the | — th frame, the m — th recovered symbol is obtained by applying a decision threshold to 
the output correlator Qı m and then the parallel symbol is converted to serial using a parallel to serial converter. 
Finally, the recovered stream bits are obtained using the Demapping function by mapping +1 to 1 and -1 to 0. 


3. PHASE CODE GENERATION 
When it comes to lowering PAPR, phase sequence selection is crucial. In this section, the phase code 
sequences are discussed. 


3.1. Random phase code 

Multiple modified data sequences are generated by multiplying random sequences with the original 
signal. The random sequences are generated by using phase sequences, which contain different rotation factors. 
The original signal phase was rotated by using these factors to reduce the PAPR. The U new generated 
sequences, which is statistically independent are forwarded into IFFT operations, then the PAPR will be 
calculated for each one to select the lowest value for transmission. The u — th generated statistically 
independent phase sequences are expressed as: 


P% = [p¥, De pi Dh ale u = 1,2,...,U (8) 


The u — th generated phase sequence is multiplied by the OFDM-DCSK symbol. The ph, represents the m — 
th phase in the u — th code, which is given by: 


P= J Pin, where Øh, E [0,27), m = 0,1,...,N — 1 and u = 1,2,...,U (9) 


The different random phase sequences are multiplied with the data symbol. The obtained sequences 
are denoted as candidates. The time domain signal is generated after taking IFFT for each candidate. After that, 
the PAPR of all time domain candidates is calculated. The lower PAPR is selected to be transmitted with the 
phase sequence which is necessary to recover the signal at the receiver. This side information is considered as 
the main drawback of this technique. 


3.2. Hadamard phase code 

In this technique, instead of random multiplying sequences, Hadamard codes were selected, which 
are perfectly orthogonal. This will improve the probability of achieving a better PAPR reduction with a small 
number of sequences. The Hadamard code is a square matrix of rows that are mutually orthogonal and entries 
which are either +1 or -1. This implies that each pair of rows in a Hadamard matrix represents two perpendicular 
vectors in geometric terms, and each pair of rows has identical entries in precisely half of their columns and 
mismatched entries in the remaining columns in combinatorial terms. As a result of this description, the related 
properties apply to all columns and rows. 


= _f7l 1 
H=[1.= [5 4l (10) 
So, the Hadamard matrix is denoted as: 

H.= A yk-1 A yk-1 | ll 
2k H k-ı —Hyk-1 (11) 

hi1 hin 
H=|: : (12) 

hyi Ann 


The m — th phase in the u“"code depending on Hadamard matrix is given by: 
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Pm = hmu M = 1,2,...,N,u = 1,2,..., N (13) 


Different U orthogonal independent vectors with a length equal to N were selected from the Hadamard 
matrix to multiply them with the data signal to generate the candidate sequences. The time domain signal is 
generated after taking IFFT for each candidate. After that, the PAPR of all time domain candidates is calculated. 
The candidate with the lowest PAPR was selected to be transmitted. The value of the index (column index) is 
needed to be transmitted as side information to recover the original data block by doing the reverse process at 
the receiver. 


3.3. Riemann phase code 

In this technique, the SLM technique is achieved by using the column phases of Riemann [25], which 
are considered as a phase sequences. The Riemann matrix is obtained by neglecting the first row and first 
column of the M x M matrix, where: 


a i—1 if idividesj 
R(ij) = 14 
GD) { —1 otherwise N 
The element value of the u — th column is either -1, which causes a phase change or 1 <u <N, 

which causes a change in the amplitude value of the modulated symbol. 


3.3. Enhanced Riemann matrix 

In this suggested technique, a MATLAB program is used to reduce the size of Riemann matrix and 
find the optimal U vectors that result in a lower PAPR and to generate a new matrix containing U columns that 
result in lower computational complexity and optimal PAPR reduction. Below is the algorithm for the 
suggested technique, and also the algorithm is shown in figure 2. 
Multiply the OFDM-DCSK signal with all Riemann column vectors. 
Taking the IFFT for all the generated sequences. 
Calculate the PAPR for all time domain sequences. 
Select the first U Riemann phase codes that have the minimum PAPR values. 
Generate the new enhanced matrix of U columns. 


eo BO oR 


N Riemann 
Phase Code 
generation 


Information 
and reference 
chaotic vectors 


Calculate 
average PAPR 
for N code 


Enhanced 
Riemann 
Matrix of U size 


Select the first U Riemann 
phase code that has the 
Minimum average PAPR 


Figure 2. Enhanced Riemann phase code based SLM technique 


4. SIMULATION RESULTS AND DISCUSSIONS 

Simulation results were taken in this section to evaluate the performance of the SLM technique by 
using Walsh Hadamard, Riemann, and the Enhanced Riemann matrix and comparing it with the conventional 
SLM technique based on computer simulations. In this work, we built the MATLAB model of the system over 
two path Rayleigh fading channel, L = 2, are used with delays, t,=0 and T} = 2T, and average power gain, 
E[a?] = = and E[a3] = > the number of subcarriers N is set to 256 and 512, the spreading factor ß is set to 
64, and the number of sequences U is selected to be 4,8, and 16. The PAPR reduction performance comparison 
of the proposed techniques was studied, while setting the SNR value to 10 dB. Figure 3 and Figure 4 describe 
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the gain of the tried techniques in terms of PAPR at CCDF = 107? & 107? for N is equal to 256 and 512 
respectively. It was noticed that the proposed technique results in a PAPR gain of approximately 9.6 dB and 
10 dB at CCDF equal to 107? when the number of subcarriers is equal to 256 and 512 respectively, which 
gives better results than other techniques for both PAPR gain value and number of candidates (U). It was 
noticed that WH SLM and conventional SLM approximately have the same PAPR value reduction, but in the 
WH technique there is no need to send side information, which is necessary in the conventional SLM technique 
to recover the original signal at the receiver. Sending the & index is enough to recover the original signal at the 
receiver for WH, Riemann, and the Enhanced Riemann techniques. Finally, the results of the performance 
evaluation are listed for all the techniques that were used in our work in Tables 1 and 2 for N=256 and 512, 
respectively. The included results are obtained for a different number of N by setting the spreading factor (6) 
to 64. The Tables 1 and 2 present the comparison between the techniques in terms of PAPRg,. 


PAPRgpf =PAPR original = PAPRrew (15) 


Table 1. Comparative study of the suggested PAPR reduction techniques for PAPRgp at N=256, p=64 


PAPRgr (dB) at CCDF=107? 


PAPRg, (dB) at CCDF=107% 


U=4 U=8 U=16 U=4 U=8 U=16 

Conventional SLM 1.75 2.3 2:7: 2.35 2.975 3.45 
WH SLM 1.75 2.275 2.7 2 2.8 3.3 
Riemann SLM 8.5 8.66 8.8 9.25 9.47 9.6 
Enhanced Riemann SLM 8.8 -- -- 9.6 -- -- 


Table 2. Comparative study of the suggested PAPR reduction techniques for PAPRgp at N=512, p=64 


PAPRgș (dB) at CCDF=107? 


PAPRgp (dB) at CCDF=107° 


U=4 U=8 U=16 U=4 U=8 U=16 

Conventional SLM 1.5 2 235 2 2.6 3.05 
WH SLM 2 2.4 2.715 2.4 2.95 3.4 
Riemann SLM 9.05 9.15 9.2 9.8 9.95 10 
Enhanced Riemann SLM 9.2 -- -- 10 -- -- 


=Pr(PAPR > PAPRo) 


CCDF(PAPRo) 


= 

o 
bo 
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Figure 3. PAPR comparisons performance for all the tried techniques for different number of candidates, 
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Figure 4. PAPR comparisons performance for all the tried techniques for different number of candidates, 
N=512 and 6=64 


5. CONCLUSION 


In this paper, SLM techniques using WH, Riemann, and the Enhanced Riemann matrix were 
performed to overcome the issue of high PAPR, which affects the efficiency of the power amplifier. The PAPR 
reduction was done by multiplying the OFDM-DCSK signal with different sequences and selecting the one 
which gives a lower PAPR value before transmitting the signal over the multipath Rayleigh fading channel. 
The inverse process was used at the receiver to recover the original signal by using the & index. The simulation 
results obtained by running the MATLAB program show that the PAPR reduction value has improved as a 
result of PAPRgr. The SLM based enhanced Riemann matrix gives the best results compared to the other tried 


techniques in terms of PAPRgp value with a smaller number of candidates U, which in turn decreases the 


processing time from 3.536 to 0.884 seconds. The elapsed time for each candidate is 0.221 seconds, which was 
calculated by MATLAB program. 
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